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Abstract 

The decay channels are investigated for extracting weak CKM phase 

information. These channels are described by color-suppressed tree diagrams only and are free from 
penguin contributions. The branching ratios for these channels are found to be ~ O (10 -5 — 10~ 6 ) 
which can be measured at the currently running B factories. The method presented here may be 
well-suited to determine the CKM angle 7. 



PACS numbers: 12.15.Hh, 13.25.Hw, 11.30.Er 



After the discovery of CP violation in K systems in 1964, there was a lot of enthusiasm 
to look for CP violation in other systems as well. CP violation is found to be rather small 
in K systems, however, theoretical predictions have suggested that the B system may be an 
ideal place for detecting sufficiently large CP violating effects. In fact large CP violation in 
B system has been observed in the currently running dedicated B factories such as BaBar 
at SLAC and BELLE at KEK. Over the years, a lot of work have been done in the area 
of CP violation with many facilities collecting data on B events. In the standard model 
(SM), CP violation arises from the non-zero weak phase in the complex Cabibbo-Kobayashi- 
Maskawa (CKM) matrix which is responsible for the charged- current weak interaction and 
is characterized by the so-called unitarity triangle. It is, therefore, imperative that the three 
angles of this triangle, termed as a, (3 and 7, be measured independently to get decisive 
information on the origin of CP violation. 

Usually, these angles are extracted from CP violating rate asymmetries in B decays. The 
angle (3 (or sin 2/3) has been cleanly determined from the measurement of the time dependent 
CP asymmetry in the golden decay mode B® — > J/ipKs [1]. The angle a can be measured 
using the CP asymmetries in B® — > 7r + 7r - |2|, but due to the existence of penguin diagrams 
there are theoretical hadronic uncertainties which are very difficult to quantify. The last 
angle which is hoped to be determined cleanly is 7 = a.rg(—V u dV* b /V c dV* b ). There have been 
many attempts, suggestions and discussions to measure this angle as cleanly as possible. 
The Gronau, London and Wyler (GLW) [2I method, which is usually termed as the golden 
method to extract 7, employs the B — > DK modes to extract this angle. In the GLW method 
the interference of two amplitudes (b — > cus and b — > ucs) is being used which endows us 
to determine the angle 7 by measuring the decay rates for B — > D°K, B — > D°K and 
B — > Dq P K and their corresponding CP conjugate modes. However, it is experimentally 
difficult to measure the mode B~ — > D°K~. As pointed out in earlier studies, the reason 
is that the final D° meson should be identified using D° — > K + tt~, but it is difficult to 
distinguish it from the doubly Cabibbo suppressed D° — > K + n~ . 

There exist many studies in the literature Q] to help overcome the difficulty and to 
provide improved ways to determine the angle 7. It is highly desirable to have independent 
measurement of the angle 7 (or otherwise the angle a), at least to the precision of the angle f3 
as of today, to understand better the CKM mechanism of CP violation under the framework 
of the SM. But so far we have not been able to succeed in this effort. Given the various 
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methods, and wide range of options available, the measurement of the angle 7 seems to be a 
better option. This is currently being done and will also be taken up in the second generation 
experiments. There is also another parameter, namely, 2/5 + 7 which is also discussed in the 
literature |5j to be measured. Since (3 is well measured by now, therefore, the measurement 
of 2/3 + 7 will be very much useful for the clean determination of 7. It should be noted 
here that we should measure the angle 7 in all possible way (and as cleanly as possible) to 
independently verify the measurements, improve the statistics and to help resolve discrete 
ambiguities. To this end, we intend to present here another important and simple way to 
extract the weak phase 7/(2/? + 7) from the decay modes 

In this paper, we consider the color-suppressed decay modes to extract 

;he CKM phase information. Several studies have been carried out using these decay modes 



for the extraction of the angle 7. We, in this investigation, present another alterna- 
tive method, which is very clean and simple to extract the weak phase tan 2 7 (tan 2 {2(3 + 7)). 
Since, as emphasized before, one needs to have as many clean methods as possible to im- 
prove the sensitivity and to resolve the discrete ambiguities the present method will be very 
much helpful in that direction. 

First let us consider the decay channels B° — ► D°K*° and B° — ► D°K*°. It has been 
shown in Ref. [3] that the CKM angle 7 can be determined by measuring the following six 
decay rates: B° -> D°K*° and B° -> D°K*°, B -> D CP K*° (where D CP = (D° + D°)/V2, 
is the CP even eigenstate of neutral D meson) and the corresponding conjugate processes. 
The D°(D°) meson is considered to decay subsequently to the flavor state K + 7t~ for which 
the ratio of the two amplitudes is found to be very tiny i.e., rp = 1^4(5° — ► K + 7v~)/A(D° — ► 
K + ix~)\ = 0.06 ± 0.003 3. Here we will show that if we consider the the decay of D 
meson to non-CP final state i.e., K* + K~ for which td ~ C(r), then it is possible to extract 
the CKM angle 7 by measuring only the four decay rates. This method is very promising 
because the experimental branching ratio for the process B° — > D°K*° has already been 
known with value Br(B° -> D°K*° = (5.3 ±0.8) x 10~ 5 and for the (B° -> D°K*°) process 
we have the upper limit as Br(B° — > D°K*°) < 1.8 x 10^ 5 [9(. The advantage of using 
the non-CP eigenstate has been discussed in 101 ]. in connection with the charged B decays 
B^ — > K ± D Q (D°), which renders the corresponding interfering amplitudes to be of same 
order. 
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Now let us denote the amplitudes for these processes as 

A B = A(B° -> D°K*°) , A B = A(B° -> D°K*°) 



(1) 



and their ratio as 
A B 



A 



r B e 



with r B 



B 



A 



B 



A 



B 



and arg (A B /A B ) = 5 B - 7 



(2) 



where 5 B and (—7) are the relative strong and weak phases between the two amplitudes. 
The ratio of the corresponding CP conjugate processes are obtained by changing the sign of 
the weak phase 7. One can then obtain a rough estimate of r B from dimensional analysis, 
i.e., 



VubV* 



VoV* 



0.4 . 



(3) 



Now we consider that both D° and D° will decay into the common non-CP final state 
(K* + K~). Denoting the D° decay amplitudes as 



A D = A(D° -> K* + K- 
one can write their ratio 



A 



D 



A(D° -> K* + K~ 



(4) 
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with td 



A 



D 



A 



D 



(5) 



where 5d is the relative strong phase between them. It is interesting to note that the 
parameters and 5p have recently been measured by CLEO collaboration JjJ, with values 
r D = 0.52 ± 0.05 ± 0.04 and 5 D = 332° ± 8° ± 11°. 

With these definitions the four amplitudes are given as 



MB° d - 


-> (K* + K~) D K*°) 


= \A B A D \ 


1 + r B r D e i{SB+SD -^ 
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MB° d - 


-> {K*-K + ) D K*°) 


= \A B A D \e i5 ° 


r D + r B e i{6B - s °- 


7)" 


MB° d - 


-> (K*~K + ) D K*°) 


= \A B A D \ 


1 + 


r B r D e i{&B+&D+l) 




MB° d - 


-> (K* + K-) D K*°) 


= \A B A D \e i5 ° 


r D + r B e i{SB - &D+l) 



(6) 



From these amplitudes one can obtain the four observables (Ri, • ■ ■ , R4), with the definition 



A t (B d ^ (K^K^d tt*°)/(A B A D ) 



(7) 
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We can now write R± = 1 + r 2 B r 2 D + 2tbTd cos(<5b + 5d — l) and similarly R 2 , -R3 and R4. 
Thus, we get four observables and three unknowns, namely, re, 5b and 7. Hence 7 can in 
principle be determined from these four observables. 

Assuming that the amplitudes \Ab\ and \Ad\ are known (so also r#, which is expected 
to be ~ 0(0.4)), we obtain an analytical expression for 7 as 

, 2 Cgl ~ -^3) 2 ~ (#2 - -R4) 2 / R s 

lan 7 [P 2 + P 4 - 2(r| + r 2 )] 2 - [Jfc + P 3 - 2(1 + r 2 r 2 ,)] 2 " W 

Thus the measurement of the four observables Ri,... t 4 can be used to extract cleanly the 
CKM angle 7. 

Next, we consider the decay channels B° d -> D*°K*°, D*°K*° and -> D*°K*°, D*°K*° 
with two vector mesons in the final state. Considering the decay of a B meson into two 
vector mesons V\ and V 2 , which subsequently decay into pseudoscalar mesons i.e., V\ — > P\P[ 
and Vi — > P2P21 one can write the normalized differential angular distribution as 121 ]. 

1 ^ J r 2/) 2/i ^2 . 2 /I • 2 n 2/ 

^3 7T3 7T~77 = ^T~F ^ L ^ cos ^1 cos ^2 + — sin 01 sin 9 2 cos ^ 

1 a cos &i a cos y 2 dip ottl I 2 

H — - sin 2 6 1 ! sin 2 9 2 sin 2 -0 H %= sin 26*i sin 26 l 2 cos ib 

2 2v/2 

h — sin 2^i sin 2^2 sin^ — — sin 2 Q\ sin 2 6*2 sin 2^ > , (9) 

2y 2 2 j 

where 9\ (62) is the angle between the three-momentum of Pi (P 2 ) in the V\ (V2) rest frame 
and the three-momentum of V\ (V 2 ) in the B rest frame, and ip is the angle between the 
normals to the planes defined by P\P[ and P 2 P 2 , in the B rest frame. The coefficients Lj can 
be expressed in terms of three independent amplitudes, Aq, A\\ and A±, which correspond 
to the different polarization states of the vector mesons as 



In = 




U 


= Re[A||A*; 


L 2 = 


l^iil 2 , 




= MA ± A* 


L 3 = 


\A±\ 2 , 


u 


= Jm[A±Al 



(10) 

In the above Aq, A\\, and A± are complex amplitudes of the three helicity states in the 
transversity basis. These observables can be efficiently extracted from the angular distribu- 



tion (jnj) using the appropriate weight functions as discussed in Ref. [131 ] . 
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The decay mode B — > V1V2 can also be described in the helicity basis, where the amplitude 



for the helicity matrix element can be parameterized as 
H x = ( Vi ( A) V2 ( A) I H e ff I B°) 



^(A)eL(A) 



hi n p 

ag ^ + p/y + e^p la p p 



(11) 



where p is the B meson momentum, A = 0, ±1 are the helicity of both the vector mesons and 
rrii, pi and £j (i = 1, 2) denote their masses, momenta and polarization vectors respectively. 
Furthermore, the three invariant amplitudes a, b, and c are related to the helicity amplitudes 
by 



H ±1 =a±cVx 2 - 1 , H = -ax - b{x 2 - 1) , (12) 

where x = (p\ ■p2)jvn\m2 = (m 2 B — m\ — mfy j {2m\m<i) . 

The corresponding decay rate using the helicity basis amplitudes can be given as 

^H \' + \H +1 \ 2 + \hA , (13) 



87rm^ 



where p cm is the magnitude center-of-mass momentum of the outgoing vector particles. 

The amplitudes in transversity and helicity basis are related to each other through the 
following relations 

A± = g+ '^ g " ' A > = A ° = H °- (14) 

The corresponding helicity amplitudes H\ for the complex conjugate decay process B — ► 
V1V2 have the same decomposition with a — > a, b — > b and c — > — c. The amplitudes a, b and 
c can be obtained from a, 6 and c by changing the sign of the weak phases. 

In order to study the feasibility of this method, first we would like to estimate the 
branching ratios of the above mentioned decay modes. Only, the experimental upper 
limits for these modes are known so far i.e., Br(B° — > D*°K*) < 6.9 x 10~ 5 and 
Br(B° — > D*°K*) < 4.0 x 10~ 5 [9]. We expect that these modes will be well measured 
soon in the currently running asymmetric B factories or in the upcoming LHCb experiment. 

In the SM, these decays proceed through color suppressed tree diagrams only and are free 
from penguin contributions. The decay B° — ► D*°K*° arises from the quark level transition 
b — ► ucs and the process B° — > D*°K*° arises from b — > cus. To evaluate the hadronic 
matrix element (Oj) = (D* ), the factorization approximation has been used. 
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Thus, in this approach, we obtain the factorized amplitude for the -Bo — > D*°K*° modes as 
H = ^A:a 2 (^ ( £l ,p 1 )|(^)^|^(p)>(^ ( £2 ,p 2 )|(«c) y _ j4 |0> 



(m B o+m K *o)A 1 (m D ,o)(el ■ e* 2 ) 
{e 1 ■ p)(e 2 • p) - «t ■ ^e fiVaf} e 2 ,1 e 1 p p? 



(15) 



(m B o + m K *o) {fn B o + m^*o) 

where / D *o is the decay constant of the vector meson D*° and A* = V* b V cs . Furthermore, 
Af K * (m^.o), A 2 K (m 2 D , ) and V B \m 2 D *o) are the form factors involved in the transition 
B° — > iT* . The coefficients a 2 is given by a 2 = C 2 + Ci/Nq, with iVc? as the number of 
colors. Thus, in this way, we can have the invariant amplitudes a, b and c (in the unit of 
G F /V2) as 

a = ia 2 X* u f D *om D *o(m B o + m K *o)A B K *(m 2 D * ) , 

u ■ \* f 2m D *om K *o BK , 2 

b = -ia 2 \ u f D *orn D *o- ■ -A 2 (m D * ) , 

(m B o + m K *o) 

■ \* f 2m D *om K *o wBK *, 2 v , . 

c = -ia 2 \ u f D *om D *v- ■ -V [m DM ) . (16) 

(m B o + m K *o) 

Substituting the values of the effective coefficient a 2 = 0.23, the Wolfenstein parameters 
A = 0.801, A = 0.2265, p = 0.189 and fj = 0.358 from Q, the decay constant f D .o = 
240 MeV, the particle masses and lifetimes from [9J] and the form factors Af K * (ra|,*o) = 
0.32,v4f^*(m^ ) = 0.31 and V BK *(m 2 D , ) = 0.52 from [16], we obtain the the branching 
ratio for the B° -> D*°K*° as 

£r(£ -> D*°K*°) = 3.87 x 10~ 6 . (17) 

Similarly, one can obtain the transition amplitude for the B° — > D*°K*° process, which is 
analogous to (IT51) with the replacement of A* by A c = V c bV* s and hence the corresponding 
branching ratio as 

Br(B° -> D*°K*°) = 2.3 x 10" 5 . (18) 

Since the branching ratios of the above two processes are very much within the reach of the 
present experiments, we expect that these processes will be measured soon by the currently 
running B factories and one will have a plenty of such events in the upcoming LHCb factory. 

Now, we consider the extraction of {2(3 + 7) from the modes Since it 

is possible to obtain the different helicity contributions by performing an angular analysis 



13l.ll7j]. from now onward we will concentrate on the longitudinal (i.e., Aq) component, which 
is the dominant one. The K s n° mode of fallows the B° -> D*°K*° and B° -> D*°K*° 
amplitudes to interfere with each other. As discussed earlier, the decay amplitude for the 
mode B® — > D*°K*° arises from b —> ucs and carries the weak phase e* 7 while B d — > D*°K*° 
arises from the quark transition b — > cms and carries no weak phase. The amplitudes also 
carry strong phases e lSl and e lS2 . Thus, we can write the longitudinal components of the 
decay amplitudes as 



Mf) 


= Amp(B d ° - 


+ /)o 


= M^V* 1 , 


Mf) 


= Amp(B° d - 


+ /)o 


= M 2 e i& , 


Mf) 


= Amp(B° d - 


+ /)o 


= M ie - i7 e i<51 


Mf) 


= Amp(5» - 


+ /)o 


= M 2 e^ 2 . 



(19) 

Since the final state / = D*°K*°is accessible to B° and B°, inserting the time evolution 
of the observables A (t) as in [18| . one arrives at the usual expression for the longitudinal 



component of the time dependent decay widths [19J as 



T (B°(t) = (\A Q (f)\ 2 + \A (f)\ 2 ) + (\A (f)\ 2 - \A (f)\ 2 ) cos Ami 

F (B°(t) - /) = £p / (|A (/)| 2 + |A (/)| 2 ) - [|A (/)| 2 - |A (/)| 2 ) cos Ami 

+ 21m ^A)(/)%(/)J sin Ami}, (20) 

where q/p = exp(— 2i[3) is the B° — B° mixing parameter and T and Am denote the average 
width and the mass difference of the heavy and light B mesons and we have neglected the 
small width difference Ar between them. 

Thus, the time dependent measurement of the longitudinal component of B°(t) — > / 
decay rates allow one to obtain the following observables : 

\A (f)\ 2 + \A (f)\ 2 , \A (f)\ 2 -\A (f)\ 2 , and Im[^ (/)M (/)], (21) 

i.e., the longitudinal components of CP averaged branching ratio, the direct CP violation 
and the mixing induced CP violation parameters. 
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Similarly one can obtain the time dependent decay rates for the final state / i.e., 
T (B°(t) - /) from by T {B°{t) -> /) by replacing A (f) by A (f) and A {f) by cor- 
responding CP conjugate A (f). T (B°(t) — > /) can be obtained from r (-B°(t) — > f) with 
similar substitution. 

Now substituting the decay amplitudes as defined in Eq. (TH?j) in (EOl) we get the decay 
rates as 

r„(B°(i)-*/) = "-j- {( Mf + Ml) + (Mf -Aff) cos Amt 

A, 



- 2M X M 2 sin(<5 - 0) sin Amt 
e -r< r 

r (5 (t) /) = — j (M 2 + M 2 2 ) - (Ml - Ml) cos Amt 



+ 2MiM 2 sin(5 + <p) sin Amt 



ro(5 u (t) -»• /) = — <j (Mf + M|) + (M{ - M|) cos Amt 



- 2M!M 2 sin(5 + 0) sin Amtj , 

r (S°(t) -»• /) = — j (M 2 + Ml) - {Ml - Ml) cos Amt 

+ 2MiM 2 sin(<5 - 0) sin Amtj , (22) 

where 5 = <5 2 — 5\ is the strong phase difference between the longitudinal components of the 
two amplitudes B° — > f and B° — > f and = 2/9 + 7. Thus through the measurements 
of the time dependent rates, it is possible to measure the amplitudes Mi and M 2 and the 
CP violating quantities S + = sin(5 + 0) and S- = sin(5 — 0). In turn these quantities will 
determine tan 2 up to a four fold ambiguity via the expression 

tan 2 0[cot 2 6} = (>S+ ~ £j= , (23) 

L J 2 - S 2 - S 2 + ± 2y / (l-5 2 )(l-5 2 ) V ; 

where one sign will give tan 2 and the other cot 2 S. 

Let us now estimate the number of reconstructed events that could be observed at the B 
factories assuming that 3 x 10 s (10 12 ) BB's are (will be) available at the e + e~ asymmetric 
B factories (hadronic B machines like LHCb). Let us first estimate the number of B° — > 
D°(D°)K* events will be available in the upcoming LHCb experiment. Assuming the 
branching ratio for the process B° -> D°K*° to be \(V ub V; s )/(V cb V* s )\ 2 xBr(B° -> D°K*°) w 



8.5 x 10" 6 , Br(D° -> K*+K~) = 3.7 x 1(T 3 [9j and 10% overall reconstruction efficiency, 
we expect to get nearly 9 x 10 3 events per year of running at LHCb. For the corresponding 
vector- vector modes, we use the longitudinal component of the branching ratio as Br (B° — > 
D*°K*°) = 0.65 x Br(B° -> D*°K*°) w 2.51 xlO" 6 , Br(D*° -> D tt ) = 62% [oj , Br(K*° -> 
Kstc ) = Br(K*° — > K7r)/3, and an overall efficiency of 10%. Thus we expect to get 
approximately 15 (5 x 10 4 ) reconstructed events at the e + e~ (hadronic) machines per year 
of running. This crude estimate indicates that this method may be well suited for the the 
extraction of the weak phase 7(2/? + 7) at LHCb. 

In this paper, we have carried out a study of the color suppressed decay modes B° — > 
D°(D°*)K°* to extract the weak phase 7(2/3 + 7). For the extraction of 7 we considered 
the decay modes B° — > D°(D°)K *, with subsequent decay of D°(D°) into non-CP state 
K* + K~ . The use of the non-CP state allows the two interfering amplitudes to be of same 
order and hence one can cleanly extract the CKM angle 7. Next we considered the processes 
B° — > D*°(D*°)K *, where the final states are admixtures of CP-even and CP-odd states. 
However, it is possible to disentangle them using the angular distributions of the final decay 
products. Now considering the longitudinal component of the time-dependent decay rates 
of these modes we have shown that </> = (2/5 + 7) can be cleanly obtained. Since these modes 
are free from penguin pollution and also the branching ratios are measurable at hadron 
factories such as the LHCb, we feel that they could be very much suited for determining the 
phase 7(2/5 + 7). 
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